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ABSTRACT 

Here we present new Keck ESI high-resolution spectroscopy and deep archival 
HST /ACS imaging for S999, an ultra-compact dwarf in the vicinity of M87, which 
was claimed to have an extremely high dynamical-to-stellar mass ratio. Our data in¬ 
crease the total integration times by a factor of 5 and 60 for spectroscopy and imaging, 
respectively. This allows us to constrain the stellar population parameters for the first 
time (simple stellar population equivalent age = 7.6+i'e Gyr; [Z/H] = — 0.95 ^q'io5 
[a/Fe] = 0.34t°;l°). Assuming a Kroupa stellar initial mass function, the stellar pop¬ 
ulation parameters and luminosity = —12.13±0.06 mag) yield a stellar mass 

of M* = 3.9to!6 x 10 6 M©, which we also find to be consistent with near-infrared data. 
Via mass modelling, with our new measurements of velocity dispersion ( a ap = 27 ± 2 
km s -1 ) and size (R e = 20.9±1.0 pc), we obtain an elevated dynamical-to-stellar mass 
ratio Md yn /M * =8.2 (with a range 5.6 ^ Md yn /M* 11.2). Furthermore, we analyse 
the surface brightness profile of S999, finding only a small excess of light in the outer 
parts with respect to the fitted Sersic profile, and a positive colour gradient. Taken 
together these observations suggest that S999 is the remnant of a much larger galaxy 
that has been tidally stripped. If so, the observed elevated mass ratio may be caused 
by mechanisms related to the stripping process: the existence of an massive central 
black hole or internal kinematics that are out of equilibrium due to the stripping event. 
Given the observed dynamical-to-stellar mass ratio we suggest that S999 is an ideal 
candidate to search for the presence of an overly massive central black hole. 

Key words: galaxies: star clusters; general - galaxies: clusters: individual: Virgo 
Cluster - galaxies: dwarf 


1 INTRODUCTION 


Ultra-compact dwarfs (UCDs) are a class of compact stellar 
system that was only discovered some 15 years ago (Hilker 
et al. 1999; Drinkwater et al. 2000. Their sizes (10 < R e < 
100 pc) and masses (10 6 < M/M© < 10 8 ) are intermedi¬ 
ate between those of globular clusters (GCs) and compact 
ellipticals (cEs). Initially, they were called ‘galaxies’ but to¬ 
day it is still debated whether UCDs are the high mass end 

* E-mail: jjanz@swin.edu.au 


of the globular cluster population (e.g. Mieske et al. 2002 
Forbes et al. 2008), or are the surviving nuclei of galax¬ 


ies that were torn apart by tidal forces (e.g. Bassino et al. 
1994 Bekki et al. 2001 Pfeffer V Baumgardt 2013). There 


is now growing evidence that both of these formation chan¬ 
nels contribute to the UCD population (e.g. Mieske et al. 
|2006| Brodie et al. |2011| Chilingarian et al. |2011[ Norris V 
Kannappan 2011 Penny et al. 2012). 

Recently Norris et al. (2014a) presented sizes, stellar 
masses and internal velocity dispersion measurements for a 
large sample of UCDs. Velocity dispersions, coupled with 
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size, provide an estimate of the dynamical mass that is rel¬ 
atively insensitive to the distribution of stellar orbits (Wolf 
et al. [2010). Early studies, based on smaller samples of 
UCDs, detected a trend for the dynamical mass-to-light ra¬ 
tio (or equivalently the dynamical-to-stellar mass ratio) to 
increase steadily with increasing stellar mass (Forbes et al. 
2008| Dabringhausen et al. 2008 Mieske et al. 12008b). Us¬ 


ing the largest sample to date of data from the literature, 
Forbes et al. (2014) showed that this trend persisted for the 
general population of UCDs even up to stellar masses of ~ 
10 9 M Q . 

Several (non-exclusive) explanations for these elevated 
mass ratios have been offered, including a non-canonical 
stellar initial mass function (IMF; both a bottom or top- 
heavy IMF would lead to an underestimation of the stellar 



per cent of the mass of the object. Although, interestingly 
despite the presence of this massive black hole M60-UCD-1 
does not display a strongly elevated mass ratio based on its 
unresolved velocity dispersion (e.g. Forbes et al. 2014 found 
a dynamical-to-stellar mass ratio of about unity). 

Forbes et al. (2014) also noted half a dozen UCDs with 
extremely high mass ratios (i.e. up to 10). Such UCDs are 
of particular interest. They may represent a rare subpopu¬ 
lation of UCDs that have been caught in the early stages 
of the stripping process, having effectively inflated sizes (or 
even outer halo structures) relative to their stellar mass. 
This may be the case for VUCD7 which reveals a core and 
halo structure in its surface brightness profile (Evstigneeva 
et al. 2008). Using the size and stellar mass of the core gives 
a dynamical-to-stellar mass ratio that is ~ 3x lower (al¬ 
though still higher than the typical UCD). We also note the 
recent discovery of a halo around the Milky Way globular 
cluster NGC 1851 (Marino et al. |2014| ). The presence of this 
halo of extra tidal light suggests that NGC 1851 was once 
the nucleus a dwarf galaxy, which was tidally stripped as it 
interacted with the Milky Way. 

However, another more mundane possibility is that 
some of the measurements in the literature for these ex¬ 
treme objects are spurious. This appears to be the case for 
M59cO which had a high mass ratio based on the measure¬ 
ment of Chilingarian & Mamon (2008) using a low resolu¬ 
tion spectrum from the Sloan Digital Sky Survey. Obtaining 
a medium resolution spectrum from the Keck telescope yield 
a new velocity dispersion measurement, which Forbes et al. 
(2014) showed gives its re-calculated mass ratio to be within 
the scatter of the UCD population as a whole. 

Hasegan et al. (2005) studied a small number of UCDs 
around M87. They calculated a dynamical mass-to-light ra¬ 
tio for half a dozen objects, finding one (called S999) to have 
the highest ratio of 9.36 =b 0.94 0 In the list of over 50 UCDs 


1 With the updated uncertainty of the velocity dispersion (see 
Section \3.3\ the uncertainty of the dynamical mass becomes 


±1.97. 


Table 1 . M87 UCDs in the field of view of the ACS observations 


Name 

Right Ascension 
( h:m:s ) 

Declination 
( d:m:s ) 

Vhel 

(km s x ) 

S999 

12:30:45.91 

±12:25:01.5 

1466 ± 5 

H44905 

12:30:57.08 

±12:23:39.8 

1563 ± 18 

S672 

12:30:54.73 

±12:21:38.3 

735 ± 106 

S887 

12:30:48.93 

±12:21:55.6 

1811 ± 106 

S928 

12:30:47.70 

±12 24 30.4 

1283 ± 5 

S5065 

12:30:50.05 

±12:24:08.9 

1578 ± 3 

S8005 

12:30:46.20 

±12:24:23.1 

1883 ± 5 

S8006 

12:30:46.65 

±12:24:22.2 

1079 ± 5 


The equinox of the coordinates is J2000. The data is taken from 
the compilation in Strader et al. (2011), following their naming 
scheme with the Hanes et al. (|2001|) extension of the Strom et al. 
( 1981) nomenclature. 


by Mieske et al. (2008b), the UCD with the highest mass-to- 
light ratio was also S999 with M/Ly = 10.2. In the Forbes 
et al. (2014) study, it has a similarly extreme dynamical-to- 
stellar mass ratio. In both the Miesk eet al.l and Forbes et al.l 
works the velocity dispersion, size and luminosity of S999 are 
based on the measurements of Ha§egan et al. (2005). 

As well as an extreme mass ratio, S999 also has a par¬ 
ticularly low stellar mass and low density for a UCD. These 
extreme properties call for a confirmation of the mass ra¬ 
tio. Here we take advantage of a very deep ACS pointing 
from the HST in two filters, to remeasure the size and pho¬ 
tometry of S999. In particular, we search for an extended 
halo in the surface brightness profile which may indicate 
ongoing stripping. We also present new high-quality Keck 
spectroscopy from which we obtain the velocity dispersion 
and stellar population properties (which in turn provide an 
improved estimate for the stellar mass). We also analyse the 
surface brightness profile of seven other UCDs around M87 
(within the deep ACS field of view) to search for extended 
halo light. 


2 OBSERVATIONS AND DATA REDUCTION 


This study is based on very deep archival HST images of 
S999 and seven more UCDs in the same field of view, and a 
newly obtained spectrum of S999 from the Keck telescope. In 
the following paragraphs the data acquisition and reduction 
are described. Here we assume that the UCDs are located 
at the same distance as M87, for which surface brightness 
fluctuations and the tip of the red giant branch give consis¬ 


tent values, i.e. 16.7 Mpc (Blake_slee et al. 2009 Bird et al. 
see also Brodie et al. 


2010 


2014) 


which means that one 


arcsecond corresponds to 1 = 80.35 pc. 


2 In the literature often a distance of 16.5 Mpc was assumed for 
the M87 UCDs. Both the inferred physical size and stellar mass 
would be reduced by 1-3%. In the mass ratio this cancels partly, 
and the derived ratio for a distance of 16.5 Mpc would be larger 
by 1.2%. 
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2.1 HST ACS Imaging 


Very deep archival HST Advanced Camera for Surveys 
(ACS) images are available in two filters (F606W and 
F814W, wide V and /, respectively) from the HST GO pro¬ 
gramme 10543 (PI Baltz). While the original purpose of the 
observations was to study microlensing towards M87, the 
images have also been used to study its globular cluster sys¬ 


tem (Waters et al. 2009 Peng et al. 2009 Madrid et ah 2009 
Webb et al. |2012 ). Similarly, we used the depth and resolu¬ 


tion of the observations to obtain surface brightness profiles 
of S999 and seven other UCDs within this field of view. Fur¬ 
thermore, we processed the overlapping ACS Virgo Cluster 


Survey pointings (ACS VCS, Cote et al. 2004) F475W ( g ) 
and F850LP (z) in the same fashion in order to enable a 
comparison of our measurements to the ACS VCS study of 
compact stellar systems (Hasegan et al. 2005). 

The pipeline processed FLT.fits images were down¬ 
loaded from the Mikulski Archive for Space Telescopes 
(MAST) archiv(0 and further reduced with the latest HST 
ACS pipeline astrodrizzle , similar to Madrid et al. (2009). 


The observations were carried out with a dither pattern al¬ 
lowing for subpixel drizzling, so that the original pixel size 
of 0705 was decreased to 07035. The drizzled and shifted 
images were combined with the IRAF task imcombine, the 
magnitudes were calibrated to the Vega system using the 
provided zeropoint^] and c orrecte d for Galactic extinction 


according to Schlegel et al. (1998) 


The field of view is filled with the light from M87. In 
order to study the light distributions of the UCDs, espe¬ 
cially at faint surface brightness levels in their outskirts, the 
light of M87 needed to be removed with high precision. This 
was achieved by first subtracting from each pixel the median 
value of those inalOlxlOl pixel box around that pixel (see 
also Madrid et al. 2009). On this preliminarily processed im¬ 
age we ran Sextractor (Bertin & Arnouts 1996) to create 
an object mask. The final model of M87 was created by a 
combination of the IRAF tasks ellipse and bmodel using 
that object mask. The final model was subtracted from the 
image before any further analysis was carried out. We found 
a remaining non-zero local background around the UCDs. 
Therefore, we subtracted the average of the mean pixel val¬ 
ues in the surrounding eight boxes of a 3 x 3 grid with a total 
width of 201 pixels centred on each UCD (see Fig. [I]). The 
uncertainty of the background was estimated by the scat¬ 
ter of these values and the mean values of 16 further boxes, 
which were placed at distances to the centre of M87 similar 
to those of the UCDs. An image of S999 and a few globular 
clusters in the surroundings are shown in Fig. [I] 


2.1.1 Surface brightness profiles 


The surface brightness profiles of S999 and seven more UCDs 
in the ACS field of view (Fig.|2|and |Al| ) were measured with 
the ellipse task of IRAF (Jedrzejewski 1987). During the 



Figure 1 . Image of S999. Left panel: F814VF-band image of S999 
with the light of M87 subtracted. The panel displays 201 X 201 
pixels or ~ 7" x 7"; at the assumed distance of 16.7 Mpc, 1" = 
80.35 pc. The grid indicates the eight surrounding boxes around 
S999 which were used to determine the local background values. 
Right panels: five globular clusters from Strader et al. (2011), 
which are close to S999 in projection, shown on the same pixel 
and grey scale for comparison. 


process fore- and background objects were masked. The po¬ 
sition angle and ellipticity of the isophotes were fixed to 
values obtained by averaging the isophotes with semi-major 
axis lengths between 0725 and 076 from a preliminary profile 
measurement (with the two quantities as free fitting param¬ 
eters) . The spatial range was chosen such that the measure¬ 
ments were unaffected by the point spread function (PSF) 
and benefited from a high signal-to-noise. The UCDs are all 
quite round (axis ratios b/a from 0.84 to 0.97). However, 
the isophotes in the outskirts can be severely biased by ex¬ 
treme pixel values due to noise in the images. The not well- 
constraint isophotes vary unreasonably to random values of 
position angle and ellipticities, and especially the latter bi¬ 
ases the surface brightness profile. For our specific question - 
looking for signs of extended envelopes in the outer parts of 
the UCDs - this was much less acceptable than the potential 
loss of information from fixing position angle and ellipticity. 
The azimuthally averaged surface brightness profiles reach 
typically a surface brightness in the /-band of ~ 24 mag 
arcsec -2 at 3 to 5 half-light radii. 

We fitted Sersic profiles to the measured surface bright¬ 
ness profiles along the semi-major axes using a nonlin¬ 
ear least-squares Levenberg-Marquardt algorithm, weight¬ 
ing the individual measurements with the inverse squared 
of the uncertainties as provided by the ellipse task. For the 
fits we excluded the innermost 0715 so as not to be affected 
by the PSF (the typical full-width at half maximum of the 
PSF for the ACS wide-field camera is ~ 07l). 


3 http://archive.stci.edu 

4 http://www.stsci.edu/hst/acs/analysis/zeropoints/zpt. 


° The difference in the attenuation as compared to Schlafly & 
Finkbeiner (2011) is smaller than 0.015 mag in all bands. 


2.1.2 Photometry 

The photometry for the eight UCDs in the ACS field of 
view was obtained in two different ways: a non-parametric 
approach and via fitted Sersic functions. 
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Figure 2. Surface brightness and colour profiles of S999. Left: from top to bottom are the surface brightness profiles and Sersic fits, the 
residuals of the profile fits, and the F606W-F814W (V — F, black) and F475W-F850LP (g — z) colour profiles (grey, with scale on right 
side). Red and blue colours represent F814W and F606W, respectively. The shaded areas in the upper two panels display the systematic 
uncertainty of the background determination, while the error bars display the measurement uncertainty as given by the ellipse task of 
IRAF. The vertical and horizontal lines in the surface brightness profile panel indicate the half-light radius and the surface brightness 
at this radius (at the assumed distance of 16.7 Mpc, — 80.35 pc). Right: Pixel values on a linear scale for F814W (top) and F606W 
(bottom) to assess the background subtraction, with the green horizontal line displaying the average pixel value at radii between 2" and 
3" after a 2.3cr clipping. 


For the non-parametric approach we followed Janz V 
Lisker ([2008) and determined (elliptical) Petrosian aper¬ 
tures (Petrosian 1976). The total flux and half-light semi¬ 
major axis length a e were then measured in an aperture with 
twice the Petrosian semi-major axis length, and corrected for 
the missing flux following the prescriptions of Graham et al. 


(2005). The parametric photometry was calculated from the 


Sersic fits to the light profiles using the formulae provided 
in Graham V Driver (2005). 


Our photometric measurements are summarised in Ta¬ 
ble]^ The deviations between non-parametric and paramet¬ 
ric photometry can be attributed to effects of the PSF, as 
the two objects with the greatest difference are those with 
the largest Sersic n, and hence the steepest profile in the 
innermost region. The uncertainty on the individual mea¬ 
surements is small, which can be seen by the consistency 
of the systematic increase of size with wavelength along 
the sequence g-V-I-z with the positive colour gradient of 
S999 (Tables [2] and Al). Furthermore, the V — I colours 
inferred from the total magnitudes from the two methods 
of measuring photometry are consistent. Our measurements 
on the ACS VCS images are very consistent with those of 
Ha§egan et al. (2005) using the same data: the typical devi- 


hts (Ha§egan et al. used King profiles) for the three UCDs 
in common is at the 5% level. In the subsequent analysis we 
consider this as the uncertainty level for the luminosities and 
radii (see also Table [3|, and adopt the average circularised 
radius from the two Sersic fits as the size R e = 20.9 pc. 


2.2 Near-infrared photometry of S999 

In addition we measured the FG-band magnitude of S999 


from archival WIRCam (Puget et al. 2004) imaging of M87 
taken as part of the Next Generation Virgo Cluster Survey- 


IR (Munoz et al. 2014). The final calibrated FG-band stacked 
image of M87 was downloaded from the CADC archive, this 
had a total integration time of 4487s and provided image 
quality of 1.31 arcseconds; sufficient for the measurement of 
a total integrated magnitude of S999 but not for any ex¬ 
amination of spatially resolved properties. The downloaded 
image was first analysed to remove the halo of M87 using 


the procedure described in Norris V Kannappan (2011). The 


ation for the magnitudes and radii derived from the profile 


total integrated FT-band magnitude of S999 was then mea¬ 
sured with SExtractor (Bertin V Arnouts 1996) yielding 
a magnitude of 17.71 ± 0.03 or Mk s — —13.40 d= 0.03 mag. 
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Table 2. Photometry of UCDs around M87 


UCD 

Non-parametric 

F606W (U) F814W (/) 

m a e ["] m a e ["} 

m 

F606W (V) 

«e ["] Me 

Sersic fit 

n m 

F814W (/) 

a e ["] Me 

n 

b/a 

PA 

[°] 

S999 

19.78 

0.29 

19.05 

0.30 

19.71 

0.26 

19.54 

1.2 

18.98 

0.27 

18.88 

1.2 

0.96 

8.8 

H44905 

21.38 

0.27 

20.65 

0.28 

21.32 

0.25 

20.94 

0.9 

20.62 

0.26 

20.28 

0.8 

0.97 

-53.1 

S672 

20.34 

0.24 

19.62 

0.24 

20.04 

0.17 

19.20 

2.1 

19.43 

0.20 

18.78 

1.8 

0.94 

-30.0 

S887 

20.60 

0.18 

19.82 

0.19 

20.12 

0.11 

18.49 

2.6 

19.45 

0.13 

18.09 

2.4 

0.96 

-54.6 

S928 

19.27 

0.34 

18.58 

0.35 

19.19 

0.30 

19.38 

1.5 

18.52 

0.31 

18.75 

1.4 

0.92 

-6.7 

S5065 

19.67 

0.21 

18.93 

0.22 

19.45 

0.17 

18.43 

1.2 

18.74 

0.18 

17.81 

1.2 

0.97 

27.9 

S8005 

19.96 

0.39 

19.25 

0.39 

19.94 

0.36 

20.31 

1.2 

19.26 

0.36 

19.62 

1.2 

0.84 

16.6 

S8006 

20.02 

0.28 

19.32 

0.28 

19.91 

0.24 

19.64 

1.3 

19.19 

0.24 

18.90 

1.3 

0.96 

16.4 


The semi-major axis length a e of the apertures containing half of the light is given in arcseconds and the surface brightness fi e of the 
corresponding isophote in mag arcsec -2 ; m is the total magnitude. The final two columns list the minor-to-major axis ratio ( b/a ) and 
position angle (PA) in degrees. At the assumed distance of 16.7 Mpc, one arcsecond corresponds to 1" = 80.35 pc. The photometry 

based on the ACS VCS pointing is summarised in Table [AT] 


2.3 Keck ESI Spectrum of S999 

A high-quality spectrum for the UCD S999 was observed 
during the night of the 2014 March 21 with the Echelle Spec¬ 
trograph and Imager (ESI, Sheinis et al. 2002) with the Keck 
II telescope with a total integration time of 3.2 hours, re¬ 
sulting in a signal-to-noise of S/N ~ 20 per Angstrom. The 
centre of the UCD was placed in the middle of the slit, which 
was aligned with the parallactic angle (PA= —73°). In to¬ 
tal seven exposures were obtained (one with an integration 
time of 600s and six with 1800s) in the echellette mode using 
a 0"75 X 20" slit, with the seeing varying between O^ffi and 
0"6. The setup results in a resolving power of R~ 5400. The 
ten Echelle orders span a wavelength range of approximately 
3900 to 11000 A. The spatial scale corresponds to 0"12 to 
0 / . / 17 per pixel from the blue to red end of the spectrum. 

From the several frames for each of the calibrations (in¬ 
ternal hat fields, bias, arc lamps) we determined master cal¬ 
ibration frames by combining the two-dimensional spectra 
using the imcombine task in IRAF. The individual science 
spectra were reduced with the pipeline MAKEE (T. Barlow). 
The software takes care of the hat helding and bias subtrac¬ 
tion. It also uses a spectrum of a star to trace the Echelle or¬ 
ders in order to extract a one-dimensional spectrum for each 
order, which is then wavelength calibrated using the known 
lines of the arc lamps (of a previous observing run using a 
similar setup). Also an error spectrum, i.e. the individual 
uncertainties of the pixels in the spectrum, is calculated. 

The hat helding does not bring the various Echelle or¬ 
ders to a common hux level. Therefore, we took the spec¬ 
trum of a velocity standard (HR3454) observed with the 
same setup and obtained the blaze function by comparing 
that spectrum to a published spectrum of the same star 


(Le Borgne et al. 2003). The Echelle orders of the cor¬ 


rected spectrum overlap. The various exposures and dif¬ 
ferent Echelle orders are finally combined to a single spec¬ 
trum in a signal-to-noise optimised way using UVES_popler 
(M.T. Murphy). The uncertainties are propagated through 
the scaling and combination processes in order to obtain a 
hnal error spectrum. 

The ESI spectrum of S999 was htted with a penalised 
pixel htting code, pPXF (Cappellari & Emsellem 2004), us¬ 
ing stellar spectra of the elodie library (Prugniel & Soubi- 


ran 200lJ version 3.1). The library contains spectra of 1388 
stars, which cover a wavelength region of 3900 A to 6800 
A at a resolution of a ~12 km s -1 (R^ 10,000), and which 
span a wide range of spectral type and metallicity. For the 
htting, the template stars were convolved to the instrumen¬ 
tal resolution (with a constant a = 21 km s -1 ), taking into 
account the templates’ resolution. The spectral range used 
for the htting was 4050 — 5550 A. The observed and best-fit 
model spectra are shown in Fig. [3] 

To estimate the uncertainties for both of the velocity 
dispersion and the Lick indices (see below), we performed 
Monte Carlo simulations: The observed spectrum was al¬ 
tered with random modifications drawn from a normal dis¬ 
tribution with a width according to the error spectrum. The 
modified spectrum was then htted in the same way, and the 
steps were repeated 1000 (for the velocity dispersion) or 500 
(for the Lick indices) times. The adopted values and their 
uncertainties are then given by the mean and standard de¬ 
viation of Gaussian fits to the distribution of the individual 
htted values. 


2.3.1 Internal kinematics 

The recession velocity and velocity dispersion were deter¬ 
mined with the best-fitting template by pPXF. We dis¬ 
carded the calcium (Ca-)triplet lines due to a lack of a 
sufficiently high resolution stellar library for this spectral 
range. For the hnal velocity dispersion and its uncertainty we 
adopted from the Monte Carlo simulation of the htting pro¬ 
cess, i.e. (Tap — 27 ± 2 km s -1 . Therefore, we obtained a re¬ 
alistic estimate for the uncertainty in the velocity dispersion 
as demonstrated in Fig. [4] Also, systematic errors from tem¬ 
plate mismatch were minimised by using a large library of 
template stars. Using velocity standard stars (a small num¬ 
ber from our observing run, and those from Ha§egan et al. 
2005) observed with the same setup, for a similar wavelength 


range as well as for the Ca-triplet, led to consistent velocity 
dispersions but larger uncertainties (for comparisons with 
the Ha§egan et al. data see Section 3.3). We also note that 


the measured velocity dispersion is well above the instru¬ 
mental resolution and Geha et al. ( |2002 ) found that it is 
feasible to reliably measure (with an accuracy of 10%) ve¬ 
locity dispersions (a ~ 18.5 km s -1 ) below the instrumental 
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Table 3. Literature values and our measurements for S999 


(1998). The measured indices are corrected to the Lick res- 


Quantity 


Literature 


This study 


Magnitude 

M g ,AB 

Mi^ab 

M Z: ab 

M v 

Mp606W(V) 

Mf814W(I) 

M Ks 


-10.791 (H05) 


-10.79 ±0.06 

-11.69 ±0.06 

-11.40 ±0.06 
-12.13 ±0.06 
-13.4 ±0.03 


Size [pc] 

Re,g 

Fe,i 

Fe,z 

Re,F606W(V) 

Re,F814W(I) 


19.9 ±0. 2 (|H 05|> 

34. 1 (Is l it 

21.9 ± 0.5~ |H05| > 


20.2 ± 1.0 

22.0 ± 1.1 
20.8 ± 1.0 
21.3 ± 1.1 


Velocity dispersion 

Cap 

Ctot 
CO 
Co 
Co 


[km s 1 ] 

23.3 ±1.3 JH05I) 


27 ±2 
26 ±2 


25.6 ± 1.4 (H05) 
23.2 ± 1.3 (M13) 
26.2 ±1.3 (MB14) 


Metallicity 

[Fe/H] 9Z 

[Fe/H] CTl 

[Fe/H] 

[Fe/H] spec 


-1.93 (H05) 
-1.38 (TT05) 
—1.40 (MI3) 


-1.27 ±0.12 


Age [Gyr] 


7.6IJ 


+ 2.0 


Mass [10 6 AL©] 


M-dyn 

22.4 ±3.0 (H05|® 

32 ±5 

Ftdyn 

21.2 (F14) 


AL* 

2.6 (F14) 

O q+0.9 
°' y -0.6 

AL* 

5.6 (M13) 

Ratios 



AL* j ' L z 

— 

1 q+0.3 
i - d -0.2 

Fi-dyn/ M* 

8.5 (F14) 

5.6...8.2...11.2 

Ft-dyn j ^AL* 

4.45 ± 0.87 (M13) 


The references in the Table are: H05 - Ha§egan et al. |2005| Sll 
- Strader et al. |2011| M13 - Mieske et al. |2013| using H05 data; 
MB14 - Mieske & Baumgardt 2014, priv. comm.; F14 - Forbes 
et al. 2014; the values were converted to our assumed distance of 
16.7 Mpc when necessary. The velocity dispersion cr tot is from the 
mass modelling, integrating over the whole object. For the mass 
ratio we list our preferred value (bold), as well as the range from 
lower/upper estimates for dynamical and stellar mass. 

resolution with a similar setup and lower signal-to-noise ra¬ 
tio. 


2.3.2 Stellar populations 

In order to derive stellar population parameters, we mea¬ 
sured Lick indices according to the definition of Trager et al. 


olution by multiplying the index measured on the observed 
spectrum Jobs with the ratio of the index measured on the 
best-fit model convolved to the observed velocity dispersion 
(Vest) and the Lick resolution 


ICorr = lobs 


rLick 

-*best 

rObs 

J best 


(i) 


The Lick resolution for each index is taken from Schiavon 


et al. (2005), and for the convolution the template resolution 
with a FWHM of 0.55 A was subtracted in quadrature. 

We interpolated the Thomas et al. (2011|) single stel¬ 
lar population (SSP) models to a fine grid (steps of 0.02 
dex in each direction; Onodera et al. |2014 ) of age, metallic¬ 
ity ([Z/H]), and a-abundance ([a/Fe]) in order to find the 
SSP-equivalent stellar population parameters of S999 via % 2 
minimisation, 


2 (-^Thomas ^Measured) 

X = --2- 


( 2 ) 


where crindex is the la uncertainty of the index derived from 
Monte Carlo simulations. 

The indices included in the analyses were H£a, Ca4227, 
G4300, H 7a , H 7f , Fe4383, C 2 4668, Mgi, Mg 2 , Mg b, Fe5270, 
Fe5335, and Fe5406. The other indices in the wavelength 
range (H* F , CN1, CN2, Ca4455, Fe4531, H^, Fe5015, 
Fe5709, Fe5782, NaD, TiOi, Ti0 2 ) were excluded since they 
are either not well calibrated in the model or significant 
outliers in the model comparisons. In the comparison of ob¬ 
served and bestfit spectrum (Fig. [3]) there are some notable 
broader deviations (e.g. around 4150 A, 4600 A, and 4800 
A). These are caused by combining the various Echelle or¬ 
ders with imperfect corrections for the blaze function, re¬ 
sulting in a slight mismatch of the spectra’s slopes in the 
overlap region of a pair of subsequent orders. It is possible 
that only the C 2 4668 and H «5 indices are affected. Omitting 
those or using a slightly different blaze correction obtained 
from a different star do not change any of our conclusions. 

We obtain: Age = 7.6±?;g Gyr; [Z/H] = -0.95tg;^; 
[a/Fe] = 0.34[Fe/H] ’= -1.27 ± 0.12 (assuming 
[Fe/H] = [Z/H] — 0.94[a/Fe], Thomas et al. [2003). The 
predicted model colours F606W (V) — F814W(I) — 0.72 
and (g — z)ab = 1.04 for an age of 8 Gyr and metallicity 
[Z/H] = —0.95 are in good agreement with the observed in¬ 
tegrated colours of F606W(V) - F8UW(I) = 0.73 ± 0.08 
and (g — z)ab — 0.90 ± 0.08 (Fig. [5]). 

The colour gradient is a hint that the stellar population 
of S999 may not be well represented by the assumption of 
a simple stellar population. In order to estimate the effect 
on the stellar mass-to-light ratio of a more complex star 
formation history we used a model grid with the combination 
of two simple stellar populations: we combined the SSPs of 
Thomas et al. ( |2Q11| with the model of a population of an 
age of 15 Gyr, [Z/H] = —2.25, and [a/H] = 0.3, with a 
mass fraction of 10% to 90% in steps of 10%. The indices in 
the combined model were calculated as luminosity weighted 


averages using the mass-to-light ratios of Maraston (2005) 


6 Obtained via a fit of a King model. 


The smallest % 2 was achieved with half of the mass in each of 
the populations. The young population was then fitted with 
an age of 4.5^'y Gyr, [Z/H] = — O.IO+q'J!?, and [a/Fe] = 
0.10l°;i°. The total stellar mass for this model is slightly 
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4200 4400 4600 4800 5000 5200 5400 
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Figure 3. Spectrum of S999 in the rest frame. The top panel shows the observed spectrum (blue line, shifted to the rest frame) and the 
bestfit template obtained with pPXF (orange line). The grey shaded regions indicate the pass bands of indices used in the determination 
of age and metallicity (omitting H lA and Mg 2 , which overlap with H lF and Mg b , respectively), while the remaining indices are indicated 
with the dashed horizontal lines. A zoom in on the H 7 line in the inset shows the quality of the fit, with the light blue shaded area 
indicating the uncertainty of the observed spectrum. The residuals are shown in the bottom panel (red line). The residual around 4530 
A is caused by a bad column in the detector. We exclude the Fe4531 index, when fitting for the stellar population parameters, since it 
coincides with some faulty pixels of the detector. 


larger than the SSP-equivalent, but not by a lot: M* = 4.1 x 
10 6 Mq and actually within the SSP stellar mass uncertainty. 


3 RESULTS AND DISCUSSION 
3.1 Light profile 

All eight analysed UCDs follow single Sersic profiles quite 
well, although two of them show deviations at a surface 
brightness around 20 to 21 mag arcsec -2 in the I band (see 
also Appendix |Aj. S999 is the only object with a small ex¬ 
cess of light in the outer parts, which is significant in both 
bands after taking measurement and background subtrac¬ 
tion uncertainties into account (Fig. [5]). However, the excess 
is at a low surface brightness level of ~ 24-25 mag arcsec -2 
in I and V, respectively (the excess contributes only a small 
fraction of the total light, i.e. ~ 2%). For comparison, the 
two Sersic components of VUCD7 (Evstigneeva et al. 2008 ) 
contribute equally to the surface brightness of gv ~ 20.6 
mag arcsec -2 at r ~ O'M. 

The size of S999 based on the Sersic fit to the surface 
brightness profile might be underestimated, since the image 
was not deconvolved, even though the innermost region was 
excluded from the fit. However, the obtained size is consis¬ 
tent with literature values (see Table |3| , with the exception 
of Strader et al. (2011), who fitted a King model with fixed 
concentration. 


At least three UCDs - S999, S887, and S5065 - have 
colour profiles that become steadily bluer towards larger 
radii over a range with a reliable colour measurement, by 
about ~ 0.15 mag. Again, the uncertainty of the background 
subtraction is the major limitation. However, the colour pro¬ 
files are consistent with the g — z profiles measured from the 


ACS VCS images. Evstigneeva et al. (2008) found similar 


positive colour gradients for the brightest UCDs in their 
sample. 

We note that the fourth cosine term profile, i.e. disky 
or boxy isophotes, does not reveal any significant features 
correlated with the surface brightness profile of S999 (not 
shown). 


3.2 Extreme mass ratio 

We determined the dynamical mass of S999 by assuming it 
is perfectly described as a Sersic function. Since the size has 
a slight wavelength dependence, we used an average value of 
those radii from the Sersic fits to surface brightness profiles 
in F606W and F814W (R e = 20.9 pc). Using the equations 
in Ciotti (1991), we deprojected the observed Sersic function 


to create a smooth three-dimensional density distribution. 
We assumed an initial mass normalisation and calculated 
the resulting projected velocity dispersion as a function of 
projected radius. From the density distribution we sampled 
10 6 stars, calculated their projected radii, and assigned each 
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Figure 4. Distribution of velocity dispersion values of S999 from 
1000 Monte Carlo simulations, simulating random artificial noise 
according to the error spectrum. We adopt the values for the 
velocity dispersion and its uncertainty as the centre and sigma 
(see box) from a Gaussian fit to the distribution (red dashed 
line). 


star a velocity along the line of sight consistent with the 
projected dispersion at that projected radius. We convolved 
each star with a Gaussian assuming a FWHM of the seeing 
(O '/6 ~ 48.6 pc). We then projected the angular dimensions 
of the slit (O'/75 x 2 / /5) to the distance of S999 and integrated 
each stellar seeing disk over the slit. We next calculated the 
observed integrated velocity dispersion within the slit, us¬ 
ing the slit integration factors as weights in the integrated 
dispersion. Finally, we scaled the mass normalisation appro¬ 
priately until the modelled integrated velocity dispersion a 


within the slit matched the observed valueljThe resulting 
total dynamical mass is Md yn = 32 d= 5 x 10 ^M©. Given the 
large number of stars used in the simulation, we found that 
the random uncertainty in the total dynamical mass is domi¬ 
nated by the uncertainty in the measured velocity dispersion 
and that sampling uncertainties are not important. 

For the conversion from luminosity to stellar mass we 
took the stellar mass-to-light ratios assuming a Kroupa IMF 
from Maraston (2005), which are valid for the Thomas et al. 
(2011) stellar population models, which we usedj^We chose 
to use the z-band, since the mass-to-light ratio depends less 
critically on the age, [Z/ H], and [a/Fe] of the stellar popu¬ 
lations at longer wavelengths. The uncertainty of the total 
luminosity in the z -band with lower S/N as compared to 
e.g. F814W is still small, and the model M/L is readily 
available. Interpolating the model stellar mass-to-light ra¬ 
tios resulted in approximately M/L z = 1.1 to 1.6 for our 


7 Our mass modelling is basically following Hilker et al. (2007). 
However, they truncated the profiles at 500 pc (H. Baumgardt, 
priv. comm.). The resulting dynamical masses are slightly smaller 
than they would be without the truncation, with the difference 
depending on the profile. For example for the Fornax UCD4 (R e = 
24.1 pc, n = 5.5, Evstigneeva et al. 2007) it is a 3% effect. 

8 Which horizontal branch morphology (blue or red) is assumed 
for the lowest metallicity models has no influence on our mass- 
to-light ratios in the z- and K s -bands. 


observed age and metallicity including the ranges of uncer¬ 
tainty. The best guess is M/L z m 1.3, which results in a 
stellar mass of M* = 3.9 x 10 6 M© (M* = 3.3 to 4.8xl0 6 M® 
for the range of mass-to-light ratios). 

Furthermore, we converted the ground-based K s mag¬ 
nitude into stellar mass, benefitting from the even longer 
wavelength. With the measured age and metallicity, Maras¬ 
ton (2005) yielded M/Lk s = 0 . 8 T 0.1 (the range comes again 
from the uncertainties in age, [Z/ H], and [a/Fe]) and a stel¬ 
lar mass of M* = 3. 8^5 x 10 6 M©. Using the Bressan et al. 


(2012 


2013) models and the procedure described in Norris 


et al. (2014b) gave a slightly larger, but consistent mass of 
M* = 4.7 ^ 0 6 x 10 6 Mq. In both cases the uncertainties are 
dominated by those of the measured stellar population age. 
Both measurements are in good agreement with the stellar 
mass derived from HST imaging (3.9 x 10 6 M©). 

The resulting dynamical-to-stellar mass ratio is ~ 8 . The 
range of possible mass ratios is rather large, due to the un¬ 
certainties in both the stellar mass-to-light ratio and the 
dynamical mass estimate (i.e. 5.6 to 11.2 with the lowest 
dynamical and highest stellar mass and vice versa). How¬ 
ever, the mass ratio remains elevated even when adopting a 
stellar mass at the high end (M* = 4.8 x 10 6 M©) of the pos¬ 
sible range and scaling the dynamical mass down to the low¬ 
est published velocity dispersion (a = 20 km s -1 ; Ha§egan 
et al. 2005 using fxcor), i.e. Md y n/M * ^ 3.7. 


3.3 Comparison to literature 


An extreme mass ratio has previously been calculated for 
S999 (Ha§egan et al. 2005 Mieske et al. 2008b Forbes et al. 


2014), all based on the ESI spectrum of Ha§egan et al. With 
improved photometry from deeper images (factor ^60 in 
total integration time), and new spectroscopy (factor ~5 in 
total integration time) allowing for additional constraints 
on the stellar population we confirm that the mass ratio is 
indeed elevated. 


Hasegan et al. (2005) measured the velocity dispersion 


of S999 in two different ways: with cross-correlation and with 
pPXF. Their two measurements are not compatible with 
each other ( a ap = 19.6 d= 1.0 km s -1 and a ap = 25.6 d= 1.4 
km s -1 , respectively), but this might be explained by an un¬ 
derestimation of the uncertainty. We note that for both their 
methods |Ha§egan et al.| used five stars, which were observed 
with the same setup. Here we made use of the much larger 
library of high resolution ELODIE stellar templates, which 
were observed with a different instrument. Repeating their 
measurement with pPXF around the magnesium triplet re¬ 
gion, we reproduce their velocity dispersion, but obtain a 
2-5 times larger uncertainty a ap = 24.5 ± 3.5 km s -1 , which 
is not accounting for systematic uncertainties. The velocity 
dispersion measurement using our own spectrum agrees well 


with that of Hasegan et al. when using pPXF (as we do). 


Given our uncertainty estimate, which is larger than that 
of |Ha§egan et al.| for their spectrum with 5 times shorter 
total integration time, our updated larger estimate for the 
uncertainty of the data with a smaller signal-to-noise ratio 
seems consistent, showing the usefulness of the Monte Carlo 
simulation for a realistic estimate. 

Our size measurements for the ACS VCS pointings are 
in good agreement with those of Hasegan et al. (2005), 


even though our approach was simpler in the sense that 
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we excluded the innermost region from the fit of the surface 
brightness profile instead of taking into account the convo¬ 
lution by the PSF as |Hasegan et al.| did. Also, the radii mea¬ 
sured on the deeper images are very consistent. All of the 
measurements follow the mild trend to larger radii at longer 
wavelengths, consistent with the colour gradient of S999. 
Strader et al. (2011) obtained a larger size using a King 
profile with fixed concentration. The Mieske et al. (2013) 
effective radius is the average of the Ha§egan et al. (2005 ) 
radii in the two filter bands, just as we use the average of 
the radii from the two deep images. 

Hasegan et al. (2005) calculated the dynamical mass 
from their King model fit to the ACS VCS surface brightness 
profile, using the central velocity dispersion, which they ob¬ 
tained from the expected velocity dispersion profile by tak¬ 
ing into account the seeing. Mieske et al. (2013) corrected the 
calculation of Hasegan et al.| to the proper King core radius 
and obtained a slightly larger dynamical mass estimate. Our 
dynamical mass for S999 is larger than both these values, 
but in agreement within the uncertainties with the corrected 
mass estimate by Mieske et al. (2013) (see Table |3|. 

The magnitudes show an excellent agreement with 
Ha§egan et al. (2005), better than our estimate of their un¬ 
certainty. Their contribution to the uncertainty of the mass 
ratio is smaller than that of any of the other parameters. 

Hasegan et al. (2005) did not calculate the mass ra¬ 
tio and quoted the dynamical mass-to-light ratio instead. 
They listed two metallicity estimates for S999, based on the 
ACS VCS g — z colour and a C — T\ colour, respectively. 
While the former is lower than our measurement, the latter 
is reasonably close. For calculating the stellar mass, Mieske 
et al. (2013) assumed an age of 13 Gyr and extrapolated 
the Maraston (2005) stellar mass-to-light ratios for S999’s 
metallicity (they used [Fe/H] = —1.4) to get M/Ly = 2.2. 
Forbes et al. polit used the stellar mass of Norris et al. 
(2014a), which used M/Ly = 1.1 (or effectively, with S999’s 
metallicity, an age of about 3.5 Gyr). While we chose to 
use the z-band instead of the U-band, the mass-to-light 
ratio in the U-band for the stellar populations we find is 
M/Lv ~ 1.7. 

As we have seen, the dynamical mass estimates for S999 
span a range of values, but agree within the uncertainties. 
The other main driver for the large range of values for the 
dynamical-to-stellar mass ratio (a factor of ~2 between the 
lowest and highest estimate) is the characterisation of the 
stellar populations. Our stellar mass-to-light ratio lies in- 
between the estimates of Mieske et al. (2013) and Forbes 
et al. (2014). Also our preferred mass ratio is in-between 
the extremes of these two studies. However, in comparison to 
those, we clearly improved the characterisation of the stellar 
populations with the ESI spectrum. Also with the better 
constrained stellar mass the dynamical-to-stellar mass ratio 
of S999 remains elevated, i.e. with the preferred value of 
Mdyn/M * = 8.2 and a range from 5.6 to 11.2. 


3.4 Elevated mass ratios in UCDs 

Our mass ratio of ~8 implies either an underestimated stel¬ 
lar mass or extra mass that contributes to the dynamical 
mass without adding light. In the following we discuss pos¬ 
sible explanations for an elevated mass ratio in S999 and in 
UCDs in general. 


3.4-1 Non-canonical initial mass function 


If the initial mass function (IMF) differs from the assumed 
one, then the inferred stellar mass changes, mostly increases. 
This is true for both top and bottom-heavy IMFs, since 
stellar remnants of massive stars add mass but no light, 
and low-mass stars contribute mass, while adding very little 
light. Both cases have been suggested to explain the elevated 


mass ratios in UCDs (e.g. Dabringhausen et al. 2008 Mieske 
et al. 


2008a). 


However, by changing the assumed IMF, the fitted SSP 
age and metallicity will also be altered. In order to esti¬ 
mate the combined effect on the mass-to-light ratio we fitted 
miles SSP models with varying IMF slope (Vazdekis et al. 


2010 with slopes from 0.3 to 3.3, solar (^-abundances). With 


increasing slope, the fitted age generally also inc reases (as 
well as x 2 ), while the metallicity decreases. The 


Vazdekis 


|et al. | models are provided for IMFs with a power-law and 
a broken power-law (with the slope varying for stars with 
M > 0.6M®, see Appendix A of Vazdekis et al. |2003| for de¬ 
tails) . Relative to the assumption of a Kroupa IMF the mass- 
to-light ratio increases by about a factor of three for the 
broken power-law IMF with the steepest slope for high-mass 
stars (and by about a factor of four for the shallowest). Us¬ 
ing a single power-law IMF, the mass-to-light increases by a 
factor of about 3 for a slope of 2.3 as compared to a Salpeter 


IMF (slope 1.3 in Vazdekis et al.|s parametrisation, with an 


increased mass-to-light ratio by about a third as compared 
to the Kroupa IMF; the shallowest single power-law IMF 
with a slope of 0.3 results in an increase by a factor 2 to 3). 
Formally, a single power-law IMF with a slope of around 3.3 
would result in a mass-to-light ratio sufficient to account for 
the dynamical mass. However, again formally, the models 
with shallower slopes should be preferred, since the models 
with slopes steeper than 2.3 let the % 2 increase dramatically 
(by a factor of more than > 24 for a slope of 3. 3)Q We note 
that such a slope is also beyond any observed IMF slope 


for the most massive galaxies (e.g. Cappellari et. al. 2012 


Spiniello et al. 2012 Conroy V van Dokkum 2012). Assum¬ 
ing that similar factors for the mass-to-light ratios would 


result from varying the IMF in the Thomas et al. (2011) 


models, which we used for our mass-to-light estimation in 
order to account for the a- enhancement, we conclude that 
none of the reasonable fits yields a mass-to-light ratio suffi¬ 
cient to result in a dynamical-to-stellar mass ratio of unity. 


We note that Frank (2014) did not find any evidence 
for a bottom-heavy IMF in S999 (nor S928), using the CO 
index. However, their uncertainties were too large to com¬ 
pletely rule out a bottom-heavy IMF. Recently, Dabring¬ 


hausen et al. (2012) claimed that UCDs have high rates of 


X-ray binaries, indicative of of an overpopulation of rem¬ 


nants of massive stars (however, Phillipps et al. 2013 find 
the opposite for Fornax UCDs). We do not find any of the 
X-ray sources of Jordan et al. (2004, with an X-ray complete¬ 
ness limit of Lx ~ 5 x 10 39 erg s -1 ) to coincide with the 
position of any of our 8 UCDs. However, this does not put 
firm constraints on the number of low-mass X-ray binaries 
and the remnants of massive stars. Thus, we conclude that 


9 This most extreme model reverts the trend of fitted age in¬ 
creasing with the IMF slope, and results in the youngest fitted 
age, possibly as well indicating that the fit is not reliable. 
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the elevated mass ratio of S999 is unlikely explained with 
a bottom-heavy IMF, but may have a contribution from an 
excessive number of stellar remnants from a top-heavy IMF. 


3.4.2 Central massive black hole 


Mieske et al. (2013) showed that on average a massive black 


hole with a mass of about 15% of the UCD stellar mass could 
explain the generally elevated mass ratios seen in UCDs, and 
is consistent with progenitor galaxies having a total mass 
of about 10 9 M©. For S999 they calculated the mass of the 
central black hole to be Mbh = 25 X 10 6 M©, several times 
the stellar mass, in order to explain their mass ratio for S999. 
This would make S999 the most black hole dominated UCD 
or galaxy known by far. 

Seth et al. ( 2014| have confirmed a massive black hole 
in M60-UCD1 (note that M60-UCDl’s mass ratio was not 
previously found to be elevated based on its unresolved ve¬ 
locity dispersion). This kind of data, however, is currently 
rare. Frank et al. ( |201 1) did not find evidence for a massive 
black hole in UCD3 (a UCD in the Fornax cluster which 
does not have an elevated mass ratio). The non-detection of 
S999 in the X-rays cannot rule out a massive black hole, but 
requires that it is non-accreting. Integral held spectroscopy 
with high spatial resolution will allow one to detect the cen¬ 
tral black hole dynamically, or to put limits on its mass. 


3.4-3 Tidal Stripping 


Forbes et al. (2014) argued with the stripping simulations of 


Pfeffer &; Baumgardt (2013) that the early stages of the tidal 


stripping process may be an explanation for the observed ex¬ 
treme mass ratios seen in half a dozen UCDs, including S999. 
The mass ratios in the stripping simulations are reduced 
when taking into account the variation of velocity disper¬ 
sion and virial coefficient (Pfeffer, priv. comm.), instead of 
assuming them to be constant like Forbes et al. (2014) did. 


While they remain elevated, this puts even stronger con¬ 
straints on the specific orbits that can cause stripping to be 
strong enough to explain the extreme mass ratios. If strip¬ 
ping is the cause of the elevated mass ratio we would expect 
the UCDs to reveal extra halo light in their surface bright¬ 
ness profiles. Our results presented here indicate that S999 
has little if any extra halo light beyond a single Sersic pro¬ 
file. Thus we conclude it is unlikely that the extreme mass 
ratio of S999 is completely due to it currently undergoing 
tidal stripping. Several other UCDs with two-component 
structures (VUCD7, UCD3, UCD5, FCC303, M60-UCD1, 
2005; Drinkwater e t al.[2 003; Chilin- 


M59cO; Ha§egan et al. 


garian & Mamon 2008 


_ Evstigneeva et al. |2008[ Chiboucas 

et al. 20lT|), possibly including two of the UCDs in our sam¬ 
ple (S928 and S8005, see Appendix[A|, are more likely candi¬ 
dates, and may suggest a continuity from nucleated galaxies 
via UCDs with envelopes to single component UCDs. 

However, stripping is still part of the explanation, if a 
central massive blackhole inherited from the more massive 
progenitor galaxy which had been stripped, is responsible 
for the elevated mass ratio. Also, the positive colour gradi¬ 
ents found for S999 and some of the other analysed UCDs 
might be expected from the stripping scenario: the assumed 
progenitors, which were nuclei of early-type dwarf galaxies, 


typically have bluer colours than the host galaxy (Lotz et al. 
Cote et al. |2006| ; Turner et al. |2012| den Brok et al. 


2004 


2011). Generally, early-type galaxies with masses of the po¬ 


tential progenitors do not show negative colour gradients 
(as seen in more massive galaxies) but often have positive 
gradients (Tortora et al. 2010), generally interpreted as age 
gradients. The remnants in the stripping simulations of Pf¬ 
effer &; Baumgardt ( |2013 ) are larger in size than the nuclei 
of the progenitor galaxies. The authors did not attribute 
this increase in size to dynamical heating of the nucleus, 
but rather to stars from the progenitor remaining in the 
remnant. As these stars are expected to be typically found 
at larger radii than the stars of the nucleus, the observed 
positive colour gradient might be expected. The typical dif¬ 
ference in colour between the nucleus and the main body of 
the galaxy (e.g. Cote et al. 2006; den Brok et al. 2011) is 
similar to the typical observed colour gradient in UCDs. 

While Evstigneeva et al. (2008) argued that the levels 
of ^-enhancements in UCDs and nuclei are not consistent 
with the stripping scenario, Paudel et al. (2010) found them 
to be similar when comparing nuclei and UCDs at similar 
local projected galaxy density. The level of (^-enhancement 
in S999 is similar to that of |Paudel et al.[ s high local density 
UCDs. 


3-4-4 Dark Matter 


Ha§egan et al. (2005) suggested (non-baryonic) dark matter 
as a possible explanation for the elevated mass ratio of S999. 


However, Murray (2009) argued that, by comparison to the 


central dark matter densities of theoretical dark matter pro¬ 
files and observed Local Group dwarf spheroidals, dark mat¬ 
ter cannot contribute large amounts of mass in UCDs. Typi¬ 
cally, the dark matter fraction within the half-light radius is 
expected to be less than 30% (see also Tollerud et al. 

Seth et al. 2014). 


2011 


In their stripping simulations, Pfeffer & Baumgardt 


(2013) assumed that the progenitor galaxy is largely 


stripped of dark matter before the stars can be affected, 
and did not include any dark matter. Baumgardt & Mieske 
(2008) showed that mass segregation in compact stellar sys¬ 


tems will expel dark matter particles from the dense central 
regions of these objects. However, the authors calculated 
that the time scale for this process is larger than a Hubble 
time for UCDs. 

We conclude that, while S999 is a good candidate to 
search for dark matter in UCDs thanks to its low stellar 
density compared to other UCDs, dark matter should not 
play a major role in explaining the elevated mass ratio. 


4 SUMMARY 

We used our measurements of size ( R e = 20.9 d= 1.0 pc) 
and velocity dispersion ( a ap = 27 =b 2 km s -1 ) based on 
high-quality data to determine the dynamical mass ( Md yn = 
32 d= 5 x 10 6 M©) of the ultra-compact dwarf (UCD) S999 
via mass modelling. Furthermore, we newly obtained stel¬ 
lar population characteristics (Age = 7.6 jjg Gyr; [Z/ H] — 
—0.951°'. 10 5 [a/Fe] = 0.341°'1°), and calculated the stellar 
mass M* = 3.9l°'g x 10 6 M© assuming a Kroupa stellar ini¬ 
tial mass function (IMF). With these new independent mea- 
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surements we derived the dynamical-to-stellar mass ratio of 
S999 Mdyn/M* = 8.2. Furthermore, we estimated realistic 
uncertainties of the measurements, and considered the pos¬ 
sible ranges of the values. We confirmed the mass ratio to be 
elevated beyond what is expected from stellar populations 
with a canonical IMF for the whole range of possible values 
5.6 < M dyn /M * < 11.2. 

This elevated mass ratio can potentially offer insights 
to the formation of S999 and other UCDs, and we dis¬ 
cussed possible, non-exclusive explanations. New measure¬ 
ments, e.g. of IMF-sensitive features of the spectrum and 
high-resolution integral field spectroscopy, are needed to ex¬ 
plore the contributions of the listed possibilities. 

We conclude that the most likely options include a cen¬ 
tral massive black hole (as recently found in M60-UCD1), 
and that the object is not in equilibrium so that the mass 
estimate based on the velocity dispersion is not adequate. 
Both might be related to S999 being the remnant nucleus 
of a tidally stripped progenitor galaxy, which could also ex¬ 
plain its positive colour gradient. Most likely a combination 
of several of the explanations is required to explain the ex¬ 
treme value of the mass ratio of S999, since none of the 
possibilities alone is likely to account for all of the excess 
dynamical mass compared to the mass in stars. 
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APPENDIX A: LIGHT PROFILES OF UCDS 
CLOSE TO M87 

In addition to S999 there are seven more UCDs in the ACS 
field of view. Here we present the surface brightness and 
colour profiles of the remaining UCDs. Most of the analysed 
UCDs follow single Sersic profiles well (Fig. |A1| ). While the 
proper subtraction of M87’s light is non-trivial, the consis¬ 
tency of the profiles and their fits in the two filter bands is 
reassuring. Two UCDs, S928 and S8005, show a clear dip of 
the observed profile in comparison to the fitted profile (seen 
as a positive residual) at a surface brightness around 20 to 
21 mag arcsec -2 in the I band. 

The colour gradients in the UCDs are mostly positive. 
Again, the uncertainty of the background subtraction is the 
major limitation. However, the colour profiles are consis¬ 
tent with the g — z profiles measured on the ACS VCS im¬ 
ages. For at least three UCDs - S999, S887, and S5065 - 
the colours become steadily consistently bluer by about ~ 
0.15 mag towards larger radii over a range where the colours 
are reliable. The fourth cosine term profile, i.e. disky or boxy 
isophotes, does not generally display significant features cor¬ 
related with the surface brightness profile - with the possible 
exception of S8005, which becomes disky in the outer parts 
(C 4 - 0.04). 

Our photometric measurements of the UCDs in the ACS 
VCS pointing are given in Table [AT] 
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Table Al. Imaging of UCDs from the ACS VCS 


UCD Non-parametric Sersic fit 
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Same as Table [2] but for the ACS VCS images in the g and z-bands. S672 is outside the field of view. 
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Figure Al. Surface brightness and colour profiles of UCDs around M87. Same as Fig. [ 2 ] The name of the UCD is given in the top left 
panel. 
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Figure Al - continued 
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Figure A1 - continued 
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Figure Al - continued 
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